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Abstract 

The Ritus' Ep eigenfunction method is extended to the case of spin-1 
charged particles in a constant electromagnetic field and used to calculate 
the one-loop neutrino self-energy in the presence of a strong magnetic field. 
From the obtained self-energy, the neutrino dispersion relation and index of 
refraction in the magnetized vacuum are determined within the field range 
m2 < eS < M^. The propagation of neutrinos in the magnetized vacuum 
is anisotropic due to the dependence of the index of refraction on the angle 
between the directions of the neutrino momentum and the external field. 
Possible cosmological implications of the results are discussed. 
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I. Introduction 



The main goal of this paper is to investigate the effects of magnetic fields 
on neutrino propagation, a topic that has recently received increasing at- 
tention. We are particularly interested in strong field effects. Its possible 
application to astrophysics, where fields of the order of 10^^ G, and even 
largei0, can be expected in supernova collapse and neutron stars, makes this 
subject worth of detailed investigation. Just to mention one of the several 
astrophysical applications of magnetic field effects in neutrino physics, one 
may recall the suggestion that the modification of the neutrino dispersion 
relation in a magnetized charged mediuml'i'l could serve to explain the high 
velocity of pulsarsi. 

Moreover, the presence of strong magnetic fields could have influenced 
the propagation of neutrinos in the early Universe and have an imprint in 
neutrino oscillations at those epochs!. The existence of primordial magnetic 
fields (of the order of 10^^ G at the electroweak scale) in the early Universe 
seems to be needed to explain the recent observations of large-scale magnetic 
fields in a number of galaxies, in galactic halos, and in clusters of galaxieJ^. 
These primordial magnetic fields could be generated through different mech- 
anisms, as fluctuations during the inflationary universei, at the GUT scalei, 
or during the electroweak phase transitionS, among others. 

Calculations of neutrino self-energies taking into account non-oerturbative 
effects of magnetic fields have been carried out in several wor , usmg 

the Schwinger methodlll. In the present paper, the Ritus' technique, which 
was originally developed for theelectron self-energy in QED in the presence 
of electromagnetic backgroundst3',ll^ is extended to the case of spin-1 charged 
particles. The Ritus' method is based on a Fourier-like transformation that 
diagonalizes in the momentum space the Green's functions of the charged 
particles in the presence of a constant magnetic field. This approachlll'S is 
particularly convenient for the strong field case, where one can constraint the 
calculation to the contribution of the lower Landau level (LLL). 

In this work we calculate the vacuum (zero temperature (T = 0), zero 
density {ji = 0)) contribution of the neutrino dispersion relation at strong 
magnetic field {ml <^ eB <^ M^, rUe is the electron mass and Mw is the 
W-boson mass). As discussed below, such a strong field can be expected to 
exist in the neutrino decoupling era. One of our main results is the existence 
of an anisotropic propagation of neutrinos in the strong magnetic field, even 
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in the absence of a medium (/i = 0). The anisotropy is due to the dependence 
of the index of refraction on the angle between the directions of the neutrino 
momentum and the external field. We also find that the terms explicitly 
depending on the mass of the charged lepton are negligible small (of order 
1/M^), while the leading term results of order thus rather significant. 

The plan of the paper is as follows. In Section II, for the sake of un- 
derstanding and completeness, we review the Ritus' method for the Green's 
function of spin-1/2 particles in the presence of a constant magnetic field. 
Then, we extend this method to the spin-1 charged particle case in the back- 
ground of a constant magnetic field (corresponding to the crossed electro- 
magnetic field (E ■ B = 0) case). In Section III, we use the results of Section 
II to calculate, in momentum space, the one-loop neutrino self-energy in the 
presence of a constant magnetic field. The neutrino dispersion relation and 
index of refraction are obtained in Section IV in the strong-field approxima- 
tion {ml <^ eB <C M^). In Section V, we make our final remarks and discuss 
a possible cosmological realization of the adopted strong-field approximation. 

II. Green's Functions at 5 7^ in the Momentum Rep- 
resentation 

The diagonalization, structure and properties of the Green's functions 
of the electron and photon in an intense magnetic field were considered ex- 
actly in external and radiative fields by Ritus in Refs. [15] and [16]. Ritus' 
formulation provides an alternative method to the Schwinger approach to ad- 
dress QFT problems on electromagnetic backgrounds. In Ritus' approach the 
transformation to momentum space of the spin-1/2 particle Green's function 
in the presence of a constant magnetic field is carried out using the Ep{x) 
function^'@, corresponding to the eigenfunctions of the spin-1/2 charged 
particles in the electromagnetic background. The Ep{x) functions plays the 
role, in the presence of magnetic fields, of the usual Fourier e*^^' functions 
in the free case. This method is very convenient for strong-field calcula- 
tions, where the LLL approximation is plausible and for finite temperature 
calculations. 

In this section we will extend the Ritus' method to the case of spin- 
1 charged particles. This extension will allow us to obtain a diagonal in 
momentum space Green's function for the spin-1 charged particle in the 
presence of a constant magnetic field. Our results are important to investigate 
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the behavior of charged W-bosons in strong magnetic fields. 

A. Electron Green's function 

For the sake of understanding, we summarize below the results obtained 
by RitusEia for the spin-1/2 case. The Green's function equation for the 
spin-1/2 particle in the presence of a constant electromagnetic field without 
radiative corrections is given by 

(7.n + m,)S(x,i/)=5(^)(x-i/) (1) 

where 

Il, = -id^-eAf, /i = 0,l,2,3 (2) 
Taking into account that 

[^(x,y),(7.n)^] =0 (3) 

it follows that S{x, y) will be diagonal in the basis spanned by the eigenfunc- 
tions of (7.11)^ 

(7.n)'vl/^(x) = V^p(x) (4) 

Since [(7.11)^ , S3] = [(7.11)^,75] = [S3, 75] = 0, one can easily find the 
eigenfunctions \E'p in the chiral representation, where S3 = ^7172 and 75 are 
both diagonal and have eigenvalues a = ±1 and x = ±1, respectively. The 
eigenfunctions are given by 

^p{x) = Ep,-^{x)Q„^ (5) 

with 00-^ the bispinors which are the eigenvectors of S3 and 75. We are 
considering the metric g^y = diag{—l, 1, 1, 1). 

In the crossed field case (E ■ B = 0), one can always select the potential 
in the Landau gauge 

A;'^' = Bxi5,2 (6) 

which corresponds to a constant magnetic field of strength B directed along 
the z direction in the rest frame of the system. The fermion eigenfunctions 
are then given by the combination 
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i?p(x) = 5^Ep.(a:)AH, (7) 

a 

where 

A(o-) = diag{6ai, 6^-1, 6^1, cr = ±1, (8) 
and the Ep^j functions are 

Ep^{x) = Ar(n)e*(P«^°+P2^'+P^^')D„(p) (9) 

In Eq. (D, N{n) = (47r |e5|)Vv^ is a normahzation factor and Dn{p) 
denotes the parabohc cyhnder functions^ with argument p = \l2 \eB\{x\ — 
H) and positive integer index 

n = n(A;,cr) = / + n = 0,1,2,..., (10) 

The integer I in Eq. ( p!0D labels the Landau levels. In a pure magnetic 
background the x dependence of the eigenfunctions Ep^^ drops away. 
The Ep functions satisfy 

7.nEp(x) = Ep{x)-iJp (11) 

where 

= (^^0, 0, -sgn (eB) ^2\eB\l,ps^ (12) 
One can easily check that these functions are both orthonormal 



J d^xEp\x)Ep{x) = {27TYf^\p-p') = 
and complete 

d^pEp{x)Ep{y) = 
Here we have used the notation Ep{x) = 



(27r)^(5fefc/5 {po -Po)S{p2 -p'2)S{p3 -p'3) 

(13) 

= {27r)^6^^\x - y) (14) 
7°E^7° and %d^p = Y,j dpodp2dp3. 
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Using the functions Ep as a new basis, we obtain, thank to the properties 



([TT| ) and ([T^), a representation of the fermion Green's function in the presence 



of a constant magnetic field which is diagonal in p 

S{p,p') = j d^xd'yEp{x)S{x,y)Ep>{y) 

= {2nrt\p-p') _^ (15) 
7.p + Trie 

The main idea of the Ritus' approach is, therefore, to use the eigenfunc- 
tions Ep{x), which correspond to the asymptotic states of the particles in the 
presence of a constant external electromagnetic field, to perform a Fourier- 
like transformation that diagonalizes the Green's functions in the momentum 
space. The advantage of the representation (|15[) is that the Green's function 
is simply given in terms of the eigenvalues (0). 

B. W-Boson Green's function 

Let us consider now the electroweak theory in the presence of a constant 
magnetic field corresponding to the potential (^. We choose the following 
gauge conditions 

Fa = d^A^ (16) 
= d^Z^ + a,M,<P, (17) 
F+ = D^'W^ + iawMw(j) (18) 
= D*^W; - lawMwcp* (19) 

with 

D^ = d^- leA'^' (20) 

In the above expressions the customary notationEl for the electroweak fields 
is used. 
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With the gauge conditions (|TB|)-(|T^ the Green's function equation for the 
W-boson takes the form 



(n^ + M^) 5;; - 2ieF^ . + (— - i)n^n, 

aw 



G,''{x,y)=S^'\x,y) (21) 



where 11^ is given in Eq. (0). 

To solve Eq. ( |2lD it is convenient to perform first a rotation in the Lorentz 
space using the transformation matrix 



1 

7! 



/ \ 

0110 

i ~i 
V V2 J 



which satisfies P ^ = P*, and the following relations 



n = n 

^^^^ a ^^a 



(22) 



(23) 



iP%)-'[tF^ P% = -i?(^3)% 
In the above expressions the notation 

n„ = (no,n+,n_,n3) 



(24) 



(25) 



n± = (Hi ± ^2) /V2 
was introduced. S3 represents the diagonal spin-one matrix 

S3 = diag{0, 1, —1, 0). 



(26) 



(27) 

After doing the rotation (p2D , the Green's function equation (^) can be 
written as 



(tf + M^) 5^ + 2eB {S3r /3 + ( i)n"n 



aw 



(28) 
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Now we can use the Feynman gaugdiHI, = 1 ^ and follow an approach 
similar to the spin- 1/2 case in order to find a diagonal in p solution of Eq. 

We start by solving the eigenvalue equation 

D%<!>l{x)=t<!>t{x) (29) 

where 



5%= {n^ + 2eBSsy 



(30) 



Because 



0, <l>^ (x) can be taken as a common eigenfunction to D 



and 5*3 . The eigenvalue equation for 5*3 is then given by 



(^3)%$ f (x) = r^$^ (x) , v = 0,±l (31) 
where rj denotes the different spin projections. From (|3l|) we can write 



(x) 



(32) 



In Eq. ( p2|) E° represents the eigenfunctions of S3 corresponding to the eigen- 
values 7] = 0, ±1, in the following way 



/ 1 \ 






and 



/0\ 




VI/ 



for r] = 0, 



/0\ 

1 




for rj = 1, and 



/0\ 



1 



for 1] = —1 
(33) 



Note that there is a degeneracy for rj = 0. 

The eigenvalue problem (pOf ) reduces now to find -Ffc^(x) from the differ- 
ential equation 



+ 2eBr] - k ) Fk^{x 



ri = 0,±l 



(34) 



From the definitions (|25| ) - (|26|) for the 11 operator in the rotated system, 
and taking into account Eqs. (|^) and (^, we can propose 



Ffc^(x) = exp {-ikoxo + ik2X2 + ik^Xs) fkr, (xi) 



(35) 
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Then, /^^ (xi) should satisfy 



j+^)/fc,(O=0 (36) 



where 

e = ^/2\eB\ (xi + k2/eB) (37) 

and 

' ^ + ^° ~ ~ ^'^'') ^^^^ 

Eq. (|36|) is the harmonic oscillator equation. Its physical solution requires 

hr,{0 ^0 for e ^ 00 (39) 

e = n + l/2, n = 0,1,2,... (40) 
From the condition (|40|) and the definition (pSD, one has 

fc^ = + A;^ + 2(n + 1 /2)eB + 2r]eB (41) 

Considering the mass shell condition k = —M^ in Eq. (|41|), we can 
write 

A;2 = (2n + I)e5-^,e5.5 + A;^ + M^ (42) 

Eq. (^) is the well-known energy-momentum relation for higher-spin charged 
particles in interaction with a constant magnetic field@. Here gs is the gyro- 
magnetic radio of the particle with spin S. For W bosons: Qs = 2. In (021) 



we can also .deutify tie so calledHa "zero-mode problem" at eB > M^. 

As known, at those magnetic fields a vacuum instability appears giving rise 
to a W-condensationii. In our calculations we restrict the magnitude of the 
magnetic field to eB < M^, thus, no tachyonic modes will be present. 

From condition (|4l|) , and taking into account that rj = 0, ±1, it follows 
that 

n=m, — rj — l, m = 0,1,2,... (43) 
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Then, from (^TJ) and (|^) we can write 

1? = -kl + kl + 2{m-l/2)eB, m = 0,1,2,... (44) 

where m are the Landau numbers of the energy spectrum of the W bosons 
in the presence of the magnetic field. 

The solution of Eq. (|36D satisfying the conditions (|39|) and ( ^Ol) is 

= 2-/2 exp(-eV4)i^n(e/2) (45) 

where Hn are Hermite polynomials. Finally, substituting this solution into 
Eq. (|35D we can write 

Fkr,{x) = Ar(n)e*(P°"''+^'2"'+P^"')D„(0 (46) 

with N{n) a normalization factor and -Dn(0 the parabolic cylinder functions. 

The functions (|3^) , together with the mass shell condition (^2|) , respec- 
tively define the wave function and energy-momentum relation for spin-1 par- 
ticles in the presence of a constant crossed electromagnetic field (E ■ B = 0). 
The study of the parallel-field case (E || B) can be found in Ref. [23]. 

Similarly to the spin- 1/2 case (Eq. (|^)), we can now form the transfor- 
mation matrix to momentum space for the W-boson Green's function, 

[Fk{x)l'= Yl F,,{x)[Q^^^J (47) 

»7=0,±1 

where the basis matrices of the Lorentz space il^'^^ are explicitly given by 

= diag{6r,fi, ^77,1, 5^ -1, 5^,o), r] = 0,±l (48) 

Notice that the only difference between the Ep{x) functions (Eq. (|^)) 
and the (x) functions (Eq. (^7])) is given through the basis of their matrix 
spaces. That is, for spin- 1/2 the transformation matrices are expanded in the 
spinorial basis A(cr), while in the spin-1 case (Eq. (|47|)) they are expanded 
in the Lorentz basis Q^. 

It can be easily shown that the (x) are orthogonal 

/.•.^.w^,w^(2.rr'(.-.o ,49) 
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with normalization factor given by A^^ = ■^47r \eB\/n\] and complete 



\d'kF^{x)Ft{y) = {2T:f5^'\x-y) 



(50) 



Using the completeness property (|50|), one can prove that the Green's 
function 



rF. ix] 



5 ^ 



-my) 



(51) 



is a solution of Eq. ( P8| ) in the Feynman gauge. 

We can use the matrix P ( ^2|) to perform a similarity transformation of 
the Green's function ( |5TD in order to represent it in the rectangular Lorentz 
space as 



GF{x,y)f,'' 



6 ^ 



d'^k 



(27r) A; 



rl/(i/) 



(52) 



where 



Vt,{x)=P%F,{x)p-^\ 



1 2no (x) 

_ 1 n+ (x) -iH- (x) 

~ 2 iH- (x) n+ (x) 

\ 2no (x) J 

and the Hq, 'H± functions are given by 



(53) 



Ho (x) = J^rn~l (x) , Hi (x) = Trn~2 [x] ± Tm {x) (54) 

In Eq. (^) Tn [x] represents the functions -Ffc^(x) evaluated at the different 
spin projections = 0,±1. That is, using the relation (^31) and evaluating 
?7 on each spin projection, we obtain the different values of n appearing as 
subindexes of Tn [x) in terms of the Landau levels m. 

From Eq. ( p2D and the orthogonality condition (^9]), the W-boson Green's 
function in a constant magnetic field can be written in the Lorentz rectan- 
gular frame as the following diagonal function of momenta 
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Figure 1: The order-g^ neutrino bubble graph. Sohd hne corresponds to the 
charged lepton Green function in a constant magnetic field and wiggly line 
corresponds to the W-boson Green function in a constant magnetic field. 



GF{k,k')/ = {27r)^f^\k-k') _J^'^ (55) 

k 

The eigenvalue k is given in Eq. (0). We conclude this section by stressing 
that the (x) matrices play for the spin-1 particle Green's function the same 
role as the Ep{x) matrices did for the spin-1/2 ones. 



III. Neutrino self-energy in the Ep — Fk representation 

It is knowni'i'lli'lli that, to lowest order, the neutrino self-energy in a 
magnetic field is given by the bubble diagram arising from the e — W loop 
(See fig. 1). 

This diagram has been calculated using a perturbative expansion in the 
magnetic field at T 7^ and /i 7^ (/x is the chemical potential of electrons) 
in Ref. [2]. Using the Schwinger proper-time method 3, which involves the 
magnetic interaction non-perturbatively, the e — W bubble has been calcu- 
lated at T 7^ and /i 7^ in Refs. [13], and in vacuum (i.e. T = and 
/i = 0) in Refs. [11,12]. When calculating the thermal contribution to the 
e — W bubble with /i 7^ or /i = oi'El, the interaction between the magnetic 
field and the charged W-bosons in the Green's function spectrum has been 
often neglected under the assumption that eB -C (in this approximation 
the magnetic field does not appear in the poles of the W-boson Green's func- 
tion). The W-boson Green's function in the above mentioned approximation, 
known in the literature as the "contact approximation," takes the form 
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^ J (^^"'^"'^^ (56) 

where 



^x,y)=exp(i^y^F^^''x,'^ (57) 

is the well known phase factor depending on the applied field@. That is, in 
the contact approximation the interaction of the magnetic field with the W- 
bosons is restricted, in the W-boson Green's function, to the phase factor 
([57|). The contact approximation must be carefully handled in vacuum (T = 
and yU = 0), since it causes severe ultraviolet divergences. To avoid those 
subtleties one can instead consider the modification of the Green's function 
of the W-boson due to the external magnetic field, on an equal foot with 
the electron, and then, only at the end of the calculation take into account 
that the W-boson mass is the largest scale in the problem. As shown below, 
this more careful approach will prove to be convenient and useful for the 
calculation of the vacuum contribution to the neutrino self-energy in the 
presence of a magnetic field. 



A. General formulation 

To calculate the neutrino self-energy in the one-loop approximation we 
start from 

S(a;, y) = ^-R^^Six, y)YGF{x, y)u (58) 

where L,R = |(1 ± 75), Gp{x, y)i, ^ is the W-boson Green's function in the 
Feynman gauge (|5^), and S{x,y) is the electron Green's function (p!5|), that 
can be expressed in the configuration space as 

^) = i 7?T4^'^ (^) -^^E,{y) (59) 

Since the neutrino is an electrically neutral particle, the transformation 
to momentum space of its self-energy can be carried out by the usual Fourier 
transform 
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Substituting with (||), (0) and (|5|) in (^) we obtain 



d^k rg^(a;)rL-(y )' 

(27r)^ I' + 



L 



(x) EJy) 



(61) 



Taking into account that the spinor matrices (§) satisfy the following 
properties 

A (±)t = A (±) , A (±) A (±) = A (±) , A (±) A (t) = 



7'A(±) = A(±)7", 7^A(±) = A(T)7' 



LA(±) = A(±)L, i?A (±) = A (±) i? (62) 

where the notation 7" = (7", 7'^) and 7-*- = (7^ 7^) was introduced, and using 
the definitions (|53|), (0), we obtain from (|6T|) 



i27r)U^'\p-p')m 



-i{p.x-p'.y) 



{2q^ ■ 7^ [lm-lA^)i:n~,,l-MA (-) + Irn-l,l-l{xK,_MA ( + )] 



+g„ ■ 7" [(/™-2/x)j;_2/i/) + ImA^)CAy)) A (-) 



+ {lm-2,l-l{x)I*^^2,l-liy) + Im,l-l{x)i:,j_,{y)) A ( + )] 



+g,e^^2-7.7' [(/.„-2,K^)C-2.(l/) - /™,/(x)/;,,(i/)) A (-) 
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+ {lm-2,l^l{x)i:^_2,l-M - Im,l-l{x)Cl-l{y)) A ( + )] } L (63) 

In Eq. (^) we used the compact notation 

laA^) =^a{x)Ei,{x) (64) 

with J^a (x) = Fkriix) and Eh{x) = Ep„{x). The subindexes a and h in (|6^ 
represent the number n, given in Eqs. (143|) and (|10]) respectively. In (|6^) the 
subindexes a and 6 were already written in terms of the Landau levels for 
the W-bosons (m) and electrons (/) with the help of Eqs. (^) and ([10|) . 

Note that, differently from the approach used in previous worksi'mll, the 
interaction between the magnetic field and the W-bosons is kept in ( |63D in 
the poles of the self-energy operator through the effective momentum k . 

Expression (|63D is the general formula for the one-loop neutrino self- 
energy in a constant magnetic field of arbitrary strength in the Ritus' ap- 
proach. 

B. Strong field approximation 

From now on, we assume that the magnetic field strength is in the range 
ml -C eB -C M^. Since in this case the gap between the electron Landau 
levels is larger than the electron mass square ( eB ^ ^l)-, it is consistent 
to use the LLL approximation for the electron (/ = 0). On the other hand, 
it is obvious that such an approximation is not valid for the W-bosons, so 
for them we are bound to maintain the sum in all Landau levels. In this 
approximation we have 

. . s —in'^ r r p-i{p-x-p' .y) 
i27c)U^'\p - p'Mp) = -f- / dVy / d'qJPd'k--^ =^ 

{q„ ■ 7" (/„_2,o(x)/;_2,,(z/) + /™,o(x)/;,o(y)) A (-) L 

+ q^e^-'^^u^' (/„_2,o(x)/;_2,o(z/) - /™,o(x)/;,o(z/)) A (-) L} (65) 

To perform the integrals in x and y in (^5|) we should take into account 
the formulas0'ii 
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d xe~ 



■ip.X T 
^ 1 



m'V[.X) 



i2ny 



(66) 



J yl\ym\ 

(67) 

where 



P± = yPi+P2, = arctan(p2M), = (68) 

ram{V,m) m'\l'\ 

Jm'v{pi.)= ^ j!(r - ^!(m^ - j)\ \^s9n{eB)p^f^''-^' (69) 

Thanks to the factor e"^?-/^, the contributions from large values of pj_ 
are exponentially suppressed in the electron LLL approximation. Thus, it is 
consistent to keep only the smallest power of pj_ in Jm'o{Pi.) so that 

Jm'oiPl.) - 5m',0 (70) 



After using Eqs. (1^), ( p7|) and ( [TOD in (pSQ , and integrating in fc, one 
obtains for the neutrino self-energy in the electron LLL approximation 



S(p) = -ig\ \eB\ j 



{qn-i\^m,25m,2 + ^m.O^m.o)^ (-) L + ^^£^''^''7^7^ (5„,25m,2 " ^m.O^m.o)^ (-) L] 

(71) 

where 

= (g„ - Pn)' + 2(m - l/2)e5 (72) 

Note that in this approximation the sum in the W-boson Landau levels 
is effectively reduced to the contribution of the two levels m = 0, 2. 
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Let us perform now the integration in the parallel momenta g,,. With 
this aim, we can use the Feynman parametrization to represent, after Wick 
rotation to Euclidean space and some variable changes, the integral in (|7ID 

as 

1 



where 

M'^ = ml+ [M^ + 2(m - l/2)eB - ml] z + pi {I - z) z (74) 

Then, performing the integrations in z and g,, and taking explicitly the 
sum in m, we arrive at 

^^^^ = (1^ " + ^2P/.e^"'^7.7l A (-) L (75) 

where 



(73) 



|e5| ^, _ ml\eB\ ^ ^M^^ , \eB\ - ml ^ ^, 4\eBf 

(76) 

The expression (p5D can be rewritten in a more convenient way using the 
following covariant form 

S(p) = aij^„ + asj^x + H + c$ L (77) 

In (^) is the four-velocity of the center of mass of the magnetized 
system (background) and 5^ = B^/ where B^ is the magnetic field in 
covariant notation B^ = \efj,up\u'^ . Notice that when a magnetic field is 
present, to form the structure of S we have to consider, in addition to the 
usual tensors p^, g^^, and e^^px, the vectors 5^ and u^. The four- velocity 
vector can be introduced in this case because the presence of a constant 
magnetic field picks up a special Lorentz frame: the rest frame (on which 
Up = (1, 0, 0, 0)) where the magnetic field is defined {upF^^^ = 0). 

In a trivial vacuum, S would depend only on the four-momentum show- 
ing no difference between longitudinal and transverse components (ai = 02) 
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in agreement with the Lorentz invariance of the system. However, when a 
nontrivial background is present, the structure of S is enriched with new 
terms related to the symmetries broken by the background. In the present 
situation, since the magnetic field introduces a special Lorentz frame, the 
four-velocity is needed to rewrite the structure in a covariant way, a situ- 
ation similar to the finite temperature casein. Moreover, because of its spe- 
cial direction in the 3-dimensional space, the magnetic field breaks one more 
symmetry: the 0(3) rotational symmetry. This new symmetry breaking is 
responsible for the appearance of the structure associated to the unit vector 
B. 

Notice that the separation between longitudinal and transverse momenta 
contributions in S (i.e. ai 7^ 02 in (|77D ), an effect normally occurring in the 
presence of a constant magnetic field, has also a covariant representation in 
terms of the basic tensors of the problem, 

A,=p'W,pW>'^^^ (78) 

^^=p''W,pW'^''^^ (79) 

where 

W^p = {u^'B^ - u'^B'') (80) 

= le.p^pW'^^ (81) 



The coefficients oi, 02, h, and c in ([77| ) are Lorentz scalars that depend 
on the parameters of the theory and the used approximation. From ([75|) - 
([76| ) one can see that in the strong-field approximation here considered their 
leading contributions in powers of are given by 

ai = T^rn-^' «2 = 0, b = aiXp, c = aiUp (82) 



{Any 



with 



i^p = p-u, Xp = P- B (83) 
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It is clear from the above equations that the longitudinal and transverse 
neutrino modes of propagation behave quite differently. This means that the 
strong magnetic field gives rise to an anisotropy in the neutrino propagation 
that is reflected in a neutrino self-energy mainly depending on the spatial 
momentum parallel to the applied field (0), (H). We point out that in the 
weak-field approximation a neutrino anisotropic propagation was found in 
Ref [11], while in Ref. [12] the splitting, characteristic in the presence of a 
magnetic field, between longitudinal and transverse momentum components, 
was absent. The last result was a consequence of the mass shell condition 
for vacuum ■y ■ p = that was imposed through out the calculation in [12]. 
Moreover, it should be notice that no linear term in B was found in [12], 
contrary to the behavior reported in [11], and to the one we found in A and 
A' from Eqs. (^2])-(^). 

For the case of neutrino propagation in a magnetized medium (/i 7^ 0), 
a self-energy structure =imikr to (0) has been reported!. However, in that 
case, the coefficients b and c are proportional to the difference between the 
electron and positron densities which are functions of the electron chemical 
potential. 

IV. Neutrino dispersion relation and index of refrac- 
tion at strong magnetic field 

Using the results ([771), (|8^) in the dispersion equation for neutrinos prop- 
agating in the external magnetic field 

det[7-p-S(p)] = (84) 
one obtains the following solution 

a;p ~ ["p^l (±1 + ai sin^ a). (85) 

In Eq. ([85|), a is the angle between the direction of the neutrino momentum 
and that of the applied magnetic field. Positive and negative signs correspond 
to neutrino and antineutrino energies respectively. 

To obtain the neutrino index of refraction n, we substitute (^5]) into the 
formula 
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to find 



n 1 — ai sin a (87) 

From Eqs. (^) and (^) it is clear that neutrinos moving with different 
directions in the magnetized space will have different dispersion relation and 
consequently, different index of refraction. That is, although the neutrinos 
are electrically neutral, the magnetic field, through quantum corrections, 
can produce anisotropic neutrino propagation. The order of the asymmetry 

2 \eB\ 

An asymmetric neutrino propagation depending on the difference between 
the number densities of electrons and positrons was previously found in a 
charged medium (yU ^ 0)i. There, the asymmetric term changes its ^ign 
when the neutrino reverses its motion. This property was suggested^! to 
be the cause of the peculiarly high velocities observed in pulsars. In our 
case, however, the asymmetric term in the dispersion relation (^) does not 
change its sign by changing a by —a. On the other hand, it is clear from 
( P7| ) that neutrinos moving along the external magnetic field {a = 0) have 
index of refraction similar to the one for the free case {n = 1), while the 
index of refraction for neutrinos moving perpendicularly to the direction of 
the magnetic field [a = vr/2), has a maximum departure from the free-case 
value. 



V. Final remarks and cosmological applications 

In this paper we have found the vacuum contribution (T = 0, /i = 0) of 
the neutrino self-energy in the strong-field regime (m^ -C eB ^ M^). The 
obtained self-energy depends only on the longitudinal neutrino momentum 
This fact is responsible of the strongly anisotropic neutrino propagation 
discussed in Sec. IV. Nevertheless, contrary to what occurs in the case with 
7^ 00, the asymmetric term in the dispersion relation maintains the sign 
when the neutrino reverses its motion. From (^) we can see that in our 
approximation the terms in E depending on the charged lepton mass mg are 
negligible small {l/M^,). 
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Our results may find applications in the physics of neutrinos in the early 
Universe. First of all, notice that the existence of strong magnetic fields in 
the early Universe seems to be a very plausible i deali, since they may be 
required to explain the observed galactic magnetic fields, i? ~ 2 x 10~^ G on 
scales of the order of 100 kpM. 

The strength of the primordial magnetic field in the neutrino decoupling 
era can be estimated from the following reasoning. Based on constraints de- 
rived from the successful nucleosynthesis prediction of primordial '^He abun- 
danceii, as well as on the neutrino mass and oscillation limits, an upper 
bound for the magnetic field produced in the early Universe prior to primor- 
dial nucleosynthesisii has been predicted. A formula for the upper bound at 
the QCD phase transition iS 

i 

Taking into account the cosmological constraint on the sum of stable neutrino 
masses ^mj,^ < 100 eV, the relation (p^) implies that at Tqcd = 200 MeV 

i 

the estimated upper limit for the primordial magnetic field is 

Bqcd < (89) 

On the other hand, taking into account the magnetic field effect of in- 
creasing n ^ p reaction rates in primordial nucleosynthesis, the magnetic 
field upper limit at the end of nucleosynthesis (T = 10^ K) is0 

Bns < WG (90) 

These values are in agreement with the equipartition principle that states 
that the magnetic energy can only be a small fraction of the Universe energy 
density. This argument leads to the relation B/T'^ ~ 2. 

Therefore, it is reasonable to assume that between the QCD phase tran- 
sition epoch and the end of nucleosynthesis a primordial magnetic field in 
the range 

ml<eB< (91) 

could have been present. 
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On the other hand, the early Universe, unhke the dense stellar medium, 
is almost charge symmetric (/x = 0), since the particle-antiparticle asymme- 
try in the Universe is believed to be at the level of 10~^° — 10~^, while in 
stellar material is of order one. Thus, to investigate neutrino propagation 
in cosmology, we would have to consider, in addition to the vacuum contri- 
bution, (me, Mvi/, ei?), of the neutrino self-energy, the = thermal 
contribution, (mg, M^r, eB, T), 

K, Mw, eB, T) = S^=° K, Mw, eB) + K, Mw, eB, T) (92) 

Notice that, when dealing with possible cosmological applications of our 
results, the vacuum neutrino self-energy calculated in Sec. Ill using the LLL 
approximation of the fermion Green function should be taken as a qualitative 
result. To understand this, we recall that a reasonable primordial magnetic 
field in the neutrino decoupling era would satisfy 

ml < eB ~ 2T^ (93) 

For such fields, the effective gap between the Landau levels is eB/T"^ ~ 
(9(1). Clearly, in this case the weak-field approximation, which would cor- 
respond to a high-temperature approximation, cannot be used because field 
and temperature are comparable. On the other hand, because the thermal 
energy is of the same order of the energy gap between Landau levels, it 
is barely enough to induce the occupation of only a few of the lower Lan- 
dau levels. Therefore, the LLL approximation, even though too radical here 
since, strictly speaking, this is not a clear-cut strong-field case, it will pro- 
vide a qualitative description of the neutrino propagation. In other words, 
the vacuum structure associated to the anisotropic propagation should still 
be present for fields eB ~ 2T^, although the coefficient c in Eq. ([77| ) may be 
quantitatively different. 

A more quantitative treatment of the neutrino propagation in a field 
eB ~ 2r^ would require numerical calculations due to the lack of a leading 
parameter. In this sense, the extension to spin-one charged particles of the 
Ritus' method developed in Section II can be very useful, since expression 
( p3| ) is a suitable representation to be used to numerically find the coefficients 
a, b and c of the general structure of the self-energy ([77|). 



22 



Therefore, one can expect that the anisotropic propagation of neutrinos 
at zero chemical potential, found in our calculations at the LLL approxima- 
tion, would be reflected in the propagation of neutrinos in the early universe 
if primordial fields satisfying the condition (|93|) were present. Even if the 
anisotropic term results small, it would account for a qualitatively new ef- 
fect. If that is the case, one can envision that the anisotropy would leave a 
footprint in a yet to be observed relic neutrino cosmic background. If such 
an effect were detected, it would provide a direct experimental proof of the 
existence of strong magnetic fields in the early Universe. 
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